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Past experience with hybrid rockets has shown that certain motor operating conditions are 
conducive to the formation of low frequency pressure oscillations, or flow instabilities, within 
the motor. Both past and present work in the hybrid propulsion community acknowledges 
deficiencies in the understanding of such behavior, though it seems probable that the answer lies 
in an interaction between the flow dynamics and the combustion heat release. 

Knowledge of the fundamental flow dynamics is essential to the basic unders tandin g of 
the overall stability problem. The present work represents a first step in this direction. A series 
of tests were conducted at NASA Marshall Space Flight Center (MSFC) on a laboratory-scale 
two dimensional water flow model of a hybrid rocket motor. Principal objectives were: (1) 
vis ualizat ion of flow and measurement of flow velocity distributions: (2) assessment of the 
importance of shear layer instabilities in driving motor pressure oscillations; (3) determination 
of the interactions between flow induced shear layers with the mainstream flow, the secondary 
(wall) throu gh-fl ow, and solid boundaries; (4) investigation of the interactions between wall flow 
oscillations and the mainstream flow pressure distribution. 

The test bed was a % scale model of an 11 inch (28cm) subscale Solid Rocket Combustion 
Simulator (SRCS) hybrid motor which has recently undergone (hot-fire) tests at MSFC. The test 
section was installed in a dual closed-loop water flow facility, originally developed in 1990 to 
investigate fluid flow through porous materials. Oxidizer, or main stream flow, was simulated 
by a water stream entering through injector slots located in the forward section of the "motor"; 
this flow could be injected either axially (along the motor axis) or radially. Fuel grain "burning” 
(mass injection only) was simulated by a secondary water stream entering the test section through 
100p sintered bronze porous plate material. The test section side walls were constructed of clear 
acrylic for flow visualization and optical velocity measurements. Flow visualization was 
accomplished through injection of very small helium bubbles into the mainstream and/or 
secondary flow streams and recording their movement with a 1000 frames/second video recorder. 

Results of these tests were presented at the 32nd A1AA/ASME/SAE/ASEE Joint 
Propulsion Conference in Orlando, FL, in July 1996, in a paper entitled Cold-Flow Study of 
Hybrid Rocket Motor Flow Dynamics". This paper (AIAA 96-2843) is included as Appendix 

1 of this report. 
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APPENDIX 1: Cold-Flow Study of Hybrid Rocket Motor 

Flow Dynamics (AIAA 96-2843) 
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Abstract 


A study is being conducted at NASA s 
Marshall Space Flight Center (MSFC) on a 
laboratory scale two-dimensional water flow 
muodelofa hybrid rocket motor. The test bed 
is a scaled model of a small hybrid motor 
which has recently undergone (hot-fire) tests 
at MSFC. The test section is installed in a 
dual closed-loop water flow facility, 
originally developed in 1990 to investigate 
fl uid flow through porous materials. 
O xidizer , or mainstream flow, is simulated by 
a water stream entering through injector slots 
located in the forward section of the "motor"; 
this flow may be injected either axially 
(along the motor axis) or radially. Fuel gram 
"burning" (mass injection only) is simulated 
by a secondary water stream entering the test 
section through lOOfi sintered bronze porous 
plate material. The test section side walls are 
constructed of clear acrylic for flow 
visualization and optical velocity 
measurements. Principal objectives include; 
( 1 ) vis ualizati on of flow and measurement of 
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flow velocity distributions; (2) identification 
of large-scale shear layer structures within the 
flow; (3) determination of the interactions 
betw ee n flow induced shear layers with the 
main stream flow, the secondary (wall) 
through-flow, and solid boundaries; (4) 
investigation of the interactions between wall 
flow o s dV h ri 1 ™* and the mainstream flow 
pressure distribution. The present paper 
describes this study and presents some 
preliminary results for (1H^)* 

Introduction 


Recent hybrid motor testing at NASA’s 
larshall Space Flight Center (MSFC) has 
town that certain motor operating conditions 
ave produced pressure traces with significant 
scillarions, or flow instabilities. Analysis of 
ata from such tests shows both high and low 
requency pressure oscillations. The high 
rcquency oscillations seem to indicate a 
pech anism driven by the first longitu dina l 
coustic mode of the combustion chamber 1 ^. 
joot pressure oscillations, on the other hand, 
yTn |0 be driven by a different mechanism 
jf set of nv-^ hg ni<ros For gaseous oxygen 
iybrid motors, r^bunixins such as chuffing, 
some boundary layer phenomenon which 
couples the combustion response to externally 
imposed pressure oscillations, and vortex 
shedding from the fuel grain face have been 
suggested as causes of these low frequency, 
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pressure oscillations' *. Other investigations 
have suggested that the flush/fill time may 
play a role in sustaining oscillations arising 
from another source 0 As a result of the 
testing done at MSFC, the question remains: 
what causes the initial oscillations? Chuffing 
has not proven to be a significant mechanism 
in the tests conducted at MSFC, primarily 
due to the nature of the oscillations and to 
post-firing grain inspection. ' The role of 
fluid dynamics with regard to initiating 
pressure oscillations has been considered, 
specifically with respect to flow fields created 
beh ind rearward facing steps as in dumps (in 
the aft cavity) and flame holders (at the 
forward end of the grain). Periodic vortex 
5 h<»dHing and downstream interactions have 
been given credit for driving pressure 
o scillations in solid propellant motors . 

One approach to resolving the conflict 
over flic driving force for these low 
frequency pressure oscillations is to better 
understand the subject the hybrid motor. A 
first step in this direction is a scries of 
jtn d ie* being conducted at MSFC on a 
laboratory-scale two dimensional water flow 
n pndgj of a hybrid rocket motor. Such water 
flow studi o have been previously performed 
elsewhere by Schadow, et. aL 1 , with regard to 
«miw phenomena associated with ramjet 
flow instabilities. 

Test objectives include; (1) 
vis ualization of flow and measurement of 
flow velocity distributions; (2) identification 
of large-scale shear layer structures within the 
flow; (3) determination of the interactions 
between flow induced shear layers with the 
main stream flow, the secondary (wall) 
through-flow, and solid boundaries; and (4) 
investigation of the interactions between wall 
flow oscillations and the mainstream flow 
pressure distribution. 


T««t Facility and Mod el Description 

The test facility was originally designed to 
study flow through porous material to 
cimtiiate solid rocket motor fuel burning in 
the Advanced Solid Rocket Motor (ASRM) 
development program at NASA’s MSFC. 
The test facility consisted of two 500 gallon 
(1892 liter) supply with one 800 gallon 

(3028 liter) return tank and a two- 
dimensional test section. Both supply tanks 
were elevated to the same height and 
regulated through water supplied from the 
return tank through a variable speed pump to 
pmintaifi equal elevation beads in each tank. 
One supply tank contained water which 
entered the test section vertically through a 
porous wall, while the other tank supplied 
water which entered the test section 
(unimpeded) horizontally. The two streams 
wall burning and mainstream gas 
flow, respectively. The water flow was 
s eeded with florcsceot dye for flow 
visualization. This concept has been a dapted 
for simulati on of the flow in a hybrid rocket 
motor. 

a . 

The configuration selected for study is 

similar to motors used in a solid rocket motor 
simulator program d uring a series of tests 
performed by Thiokol to investigate 
suppression of pressure oscillations in 
subscale hybrid rocket motors. 1 This 
configuration consisted of three main circular 
port (CP) sections of fuel butted together 
with a forward section (upstream plenum) 
which was either lined or unlined with fbel 
and an aft section (downstream plenum) 
which was imbued. The ability to inject 
oxidizer radially into the aft end was 
incorporated. Post-firing weighing of each 
of the circular port sections revealed a 
variation in the oxidizer-to-fuel (O/F) ratio 
along the motor (as expected). This aspect of 
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the hybrid was included in the model. 
Though several injectors were used, only 
radial injection patterns and axial injection 
patterns are incorporated in the model. Also 
included in Thiokol’s study of pressure 
oscill ati on suppression, as well as in the 
m od el, were inhibitor plates or rearward 
facing steps of various heights which would 
inhibit burning at the face of the fuel grain 
and change the flow pattern at the head end 
of the motor.* This study included both 
"stable" and "unstable" motors; unstable 
motors had pressure oscillations higher than 
5 % of the average motor pressure. 

The original constraints for the hybrid 
flow simulation model were set to 
accommodate minimal changes in the test 
facility while maintaining the ability to 
adequately mod*' the entire geometry of a 
hybrid motor (albeit two-dimensionally). 
This led to a m? T<m,tm test-bed length of 
approximately 36 inches (92 cm). In order to 
visualize and record phenom e na in the entire 
test section, a minimum CP flowpaih height 
(inner wall to inner wall) was set at 2 inches 
(5 cm). Based on data gathered from original 
in the ASRM model, 1/4 inch (0.64 
'em) th ick lOOp sintered bronze was selected 
to mass addition through the walls 

(top and bottom). A preliminary model of 
the test facility and test rig was used to 
derive flow losses in the system and to 
jc^rrmfia achievable Reynolds number (Re) 
levels as compared to typical Reynolds 
numbers found in actual hybrid motors. The 
results are shown in Fig- 1, plotted as Re vs. 
scaling factor. As expected, strict Reynolds 
number similarity is not achievable with this 

model. 

Constraints and requirements for the 
hybrid flow simulation model, as determined 
by the test objectives, were as follows: 


1. easily changeable configurations 

2. flaw visualization capability 

3. laser doppler velodmetry (LDV) 
capability 

4. particle displacement tracking (PUT) 
capability 

5. O/F ratio similar to actual motor 

6. geometry simulation of actua l motor 
including forward and aft steps 

7. various sections of fuel burning ’Anass 
addition (3 piece circular port and head- 
end lining) 

8. various areas of *oxidizer" injection 
(injector and aft end) 

9. oscillating secondary flow (to assure an 
" unstable" configuration) 

10. inlet and exit areas simulating actual 
motor (injector and nozzle geometry) 


NASA/MSFC Dual-Loop Water Test Rig 



Fig. 1: Facility Reynolds number study 


The model scaling rationale was as follows: 


• Initial scale was based on the actual motor 
size and the test facility limits far model 
length. 

• Model size chosen was a trade-off between 
Reynolds number scaling (small scale 
desired) and ease of visualization and 
measurement (large scale desired). 

• The scale factor and the CP diameter of 
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the actual motor fixed the model CP 
hydraulic diameter, 

• A 2.0 inch x 0.875 inch (5 cm x 2.2 cm) 

rectangular CP passage was selected to 
satisfy the hydraulic diameter criterion. The 
5 cm height was chosen for ease of 

flow measurement. 

• The upstream plenum configuration 
resulted from the constraints of hydraulic 
diameter (2.33 inches or 5.9 cm), step height 
scaling, and scaling of Reynolds number 
chan ge between the plenum and the CP 
portion. The result was a 3.12 inch x 1.86 
in ch ( 7.9 cm x 4.72 cm) rectangular passage. 

• The downstream plenum has the same 
(througbflow) dimensions as the upstream 
plenum. 

The lengths for the upstream plenum. CP 
portion, downstream plenum, and n ozzle 
ramp out sections were scaled directly from 
the Qf *t ual mo tor dimensions. The model 
nozzle angle was die same as in the actual 


motor. The injector opening area was not 
scaled from actual motor injector data due to 
l irmtatinns in the available test facility 
pressure head. Since the injector produces 
the largest flow restriction, injection area was 
determined from O/F requirements and 
pressure drop. The ratio of the required open 
area of the inhibitor plates to the CP flow 
area was determined based on the size of the 
inhibitor plates used in the Thiokol tests. 
Three inhibitor plates were made, with open 
area ratios of 20%, 36% and 64% (as was a 
fourth plate, with 100% open area). To 
provide a mean” to externally oscillate the 
secondary or wall flow, a modified systolic 
p ump was installed in the feed line for that 
loop. The pump was connected to a variable 
speed drive which allowed for forced 
frequencies from 6 Hz to in excess of 25 Hz. 
None of the results reported here utilized the 
external oscillation system, however. 
Schematics of the model and test facility are 
shown in Figs. 2 and 3, respectively. 





Fig. 2: Model Schematic 
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Fi g . 3: Test Facility Schematic 

Tnst rif flffiflation end Test Procedure 

Four test configurations were studied in 
this phase of the project. Confiptrdlion l t 
the baseline configuration, consisted of the 
model with the primary mass flow 
("oxidizer") injected axially into the upstream 
plenum and secondary mass flow ("fuel") 
injected through the upstream plenum walls 
and through all sections of the CP portion. 
The 100% open plate was used (no inhibitor), 
and no aft end (downstream plenum) 
secondary wall flow was utilized. 
Prel iminar y testing for this configuration 
Included seeding the flow with laser sensitive 
dye and recording the head and aft sections at 
both 500 and 1000 frames per second (fps)- 
The upstream plenum secondary flow was 
then seeded with helium bubbles and back¬ 
lighting was used to do flow visualization at 


1000 fps. Following this, LDV 

measurements were marie. 

«% 

Co 9 \figvf 8 tiofi 2 utilized the baseline 
configuration except the radial injector was 
incorporated. Testing included flow 

visualization with laser sensitive dye, flow 
visualization with helium bubbles, and LDV 
measurements* 

In configuration 3 the 64% open area flow 
inhib itor plate was installed in the baseline 
configuration (axial injector)* Testing with 
this configuration included flow visualization 
with helium bubbles and LDV m ea s u r e m e nt s* 
It was determined after processing the flow 
viffusli z^tinn data that the laser sensitive dye 
added no real value to the study being 
conducted* 
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Configuration 4 was the same as 
configuration 2 (radial injector) except that 

the 64% open area flow inhibitor plate was 
installed. Flow visualization with helium 
bubbles and LDV measurements were made. 

LDV measurements were made with 
a TSI Inc. fiber optic based system. Shown 
schematically in Fig. 4, the system uses a 300 
miUi'Watt (mW) Argon Ion Laser that results 
in approximately 50 mW of laser power in 
the measuring volume. The optical signal is 



amplified and processed by TSI's ColorLink 
and IFA 750 respectively. A fifteen 
m illime ter diameter probe was used, 
providing excellent measurement volume 
spatial resolution (see Fig. 4). Alu min u m 
oxide particles of 1 micron diameter wen; 
a dded to the water supply tank for LDV 
seeding. These particles are excellent for 
water because they mix easily and are 
inexpensive. (Union Carbide Corporation 
Specialty Powders Business, 1555 Main Street. 
Indianapolis, IN 46224). 

Velocity measurements were made in both 
the upstream plenum and CP portions of the 
model. The measuring volume was centered 
in the model by visual inspection, and probe 
movements were acc ompl is h ed with computer 
control of a traverse table. Excellent data 
rates were maintained by stirring up the 
particles that settled in the storage tanks. A 
hi gh speed video system, NACs HSV-1000, 
was used to generally characterize the 
upstream plenum and CP flow fields. He li u m 
bn hhfc* were injected Into the oxidizer 
flowstream and illuminated from the rear. In 
a/iH«ri«n to the LDV measurements, video 
pictures w e re made at 500»1(J00 fps. 

P«CT»ltS 

Slow motion replay of the high-speed video 
dearly shows generation of vortices (in most 
instances) from the inhibitor simulator and 
flow recirculation zones in the bead end for 
both radial and axial injectors. Shedding 
as obtained by visual count of 
vortices over a specified time interval (~ 0.4 
yco nd*. real time), are as follows: 

Configuration 1: not discemable 
Configuration 2: 35 Hz 
Configuration 3: 30 Hz 

Configuration 4: 39 Hz 


Fig. 4: LDV system schematic 
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Estimated uncertainty in the above values is ± 
20 %. for configuration 1 were too 

uncertain to report These frequencies arc in 
the expected range for shed vortex pairs with 
a Strchal number of the order 0.25 - 0.50. 

Figures 5 -12 show the LDV results for each 
configuration. In each case, approximately 
1200 points were taken (both u and v 
velocity components). Oat* wot then 
interpolated and plotted in either streamline 
foaaat or as absolute velocity contours (the 
square root of the sum of the squares of the 
velocity components). 



5; Configuration 1 LDV streamlines 



Fig. 6: Configuration 1 LDV velocity 
contours 



Fig. 7: Configuration 2 LDV streamlines 



Fig. 8: Configuration 2 LDV velocity 
contours 



fig. 9: Configuration 3 LDV streamlines 
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Fig. 10: Configuration 3 LDV velocity 
contours 


Clearly seen in the data are such features as 
the injected mainstream flow, the injected wall 
flow, and, in some instances, large-scale 
recirculation regions. 

Finally, Figures 13 and 14 show the 
jpfliignce of the 64% open inhibitor plate on 

the overall (time averaged) flow. This is done 

by taMng the difference between the data sets 
obtained for the 100% open and 64% open 
cases and plotting the results as absolute 
velocity contours. Local velocity differences 
range up to about 9 ft/sec (2.75 m/sec). For 
clarity, differences less than 1 ft/sec (03 
m/sec) are not shown. 




fig. 11: Configuration 4 LDV streamlines 


Fig. 13. Effect of inhibitor, axial injection 



pig. 12: Configuration 4 LDV velocity 
contours 


Fig. 14: Effect of inhibitor, radial injection 
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Summary and Conclusions 


An understanding of the basic flow 
dynamics of the subject hybrid model has 
been gained through this series of “sting. 
r~T-ianginp injectors (axial vs. radial) and 
inhibiting the flow between the upstream 

plenum and the CP section changes the basic 

flow structure, as evidenced by streamline and 

velocity contour plots. Numerous shear yer 
structures were identified in the test 
configurations; these structures include both 
yqluting and traveling vortices which may 
affect combustion stability- Standing vortices 
may play a role in the heal addition process as 
oxidizer enters the motor, while traveling 
vortices may be instability mechanisms m 
themselves. Finally, the flow visualization and 
LDV measurements give insight into 
determining the effects of flow inducedshear 
layers. One such case is illustrated in Fig. 7 
where the CP wall flow (and possibly the 
Hajot flout in a hot flow situation) is being 
drawn into the upstream plenu m , whereas Fig. 
11 shows that the presence of an inhibitor may 
reduce this tendency. 


When the external oscillation system 
was used in configurations 1 and 4, no 
changes in the flow characteristics could be 
detected visually. The LDV measurements, 
being time-average values, did not discern any 
changes either. Therefore, one modification to 
the will involve the addition of fifteen 

n s> fluctuating pressure transducers in various 

locations on the model. This, coupled with 

(be inclusion of a 45“ injector and the 20% 
and 36% inhibitor plates, will further extend 
tbe understanding of the model flow 
characteristics. Particle displacement tracking 
(PDT) will also be done with configurations 1- 
4 and with subsequent test configurations 
The results of all testing will be compiled and 
used to build a CFD model to predict flow 


behavior. Once the CFD model is in place, 
hypotheses linking flow dynamics to hybrid 
motor stability can be tested. 
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